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INTRODUCTION 


N'l  Design  and  evaluation  of  radio  frequency  (RF)  systems  that  must  operate  through  ionospheric 
disturbances  resulting  from  high  altitude  nuclear  detonations  requires  an  accurate  channel  model.  Such 
a model  can  be  used  to  construct  realizations  of  the  received  signal  for  use  in  digital  simulations  of 
trunsionospheric  communications  links  and  radars  or  for  use  in  hardware  channel  simulators,  j 


It  is  well  known  that  the  first  order  statistics  of  the  received  RF  signal  after  propagation  through 
strongly  scattering  ionospheric  disturbances  are  accurately  described  by  the  Rayleigh  probability 
distribution  for  the  amplitude  and  a uniform  distribution  for  the  phase.  Equivalently,  the  two 
orthogonal  components  of  the  received  electric  field  are  independent,  zero  mean  normal  random 
processes.  Given  the  second  order  statistics  of  the  received  signal,  standard  statistical  techniques  can  be 
used  to  generate  realizations  of  the  channel  impulse  response  function  from  which  the  received  signal 
can  be  constructed. 


The  second  order  statistics  of  transionospheric  RF  scintillation  are  described  by  the  two-position, 
two-frequency,  two-time  mutual  coherence  function  of  the  received  signal.  The  derivation  of  the 
mutual  coherence  function,  which  is  a solution  of  Maxwell's  equations,  requires  a model  for  the 
temporal  and  spatial  variations  of  the  electron  density  fluctuations  in  the  ionosphere.  Under  Taylor's 
frozen-ih  hypothesis,  the  electron  density  fluctuations  are  described  as  a rigid  structure  that  drifts  past 
theline-of-sight.  There  is  then  a deterministic  relationship  between  spatial  and  temporal  fluctuations  in 
the  electron  density.  This  model  is  accurate  when  the  ionization  has  formed  a thin  layer  of  striations 
/ aligned  with  the  geomagnetic  field  lines.  Before  striations  have  formed  or  when  there  are  multiple 
scattering  layers  in  the  ionosphere  with  different  relative  velocities,  a turbulent  model  may  be  more 
appropriate.  In  the  fully  turbulent  case,  the  spatial  and  temporal  fluctuations  of  the  electron  density  are 
uncorrelated.  Reality  should  lie  somewhere  between  these  two  limiting  models. 

^ This  paper  describes  a general  model  which  varies  smoothly  between  the  frozen-in  and  turbulent 
models.  The  starting  point  for  the  general  model  is  a description  of  the  electron  density  fluctuations  in 
the  ionosphere  that  includes  space-time  correlation.  Then  the  mutual  coherence  function  for  the 
general  model  is  presented.  Examples  are  given  which  illustrate  the  effects  on  the  received  signal  due 


to  the  variation  from  the  turbulent  to  frozen-in  models. 
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MUTUAL  COHERENCE  FUNCTION 

The  derivation  of  the  mutual  coherence  function  starts  with  Maxwell's  equations  from  which  the 
parabolic  wave  equation  is  derived.  The  parabolic  wave  equation  can  be  solved  to  give  the  received 
electric  field  for  a specific  electron  density  distribution  in  the  ionosphere.  However,  the  electron 
density  distribution  is  a random  process  so  the  received  electric  field  is  also  a random  process.  The 
parabolic  wave  equation  is  therefore  used  to  derive  an  equation  for  the  two-position,  two-frequency, 
two-time  mutual  coherence  function  of  the  electric  field,  r(Ar,A{i),At).  The  solution  of  the  differential 
equation  for  T,  which  is  also  a solution  of  Maxwell's  equations,  then  provides  a description  of  the 
second-order  statistics  of  the  received  electric  field.  The  Fourier  transform  of  the  mutual  coherence 
function  is  the  Generalized  Power  Spectral  Density  (GPSD)  of  the  received  signal. 


Consider  a monochromatic  spherical  wave  with  an  electric  field  E(r,(0,t)  which  is  a function  of 
position  r,  carrier  radian  frequency  to,  and  time  t.  The  wave  originates  from  a transmitter  located  at  r 
a (0,0, -zt)  and  propagates  in  free  space  in  the  positive  z direction  until  it  is  incident  on  an  irregularly 
ionized  layer  which  extends  from  0 < z < L and  is  infinite  in  the  x-y  plane.  After  emerging  from  the 
layer  at  z « L,  the  wave  propagates  in  free  space  to  a receiver  located  at  r = (0,0, zr).  This  geometry  is 
illustrated  in  Figure  1.  As  the  wave  propagates,  its  phase  varies  as  itot  - i(!c)z,  so  we  can  write 


E(r,to,t)  = U(r,to,t)  exp 


~(  z 
i <ot  - } (k(z'))  dz' 

- V *z«  > 


where 


0) 


(2) 


c is  the  speed  of  light  in  a vacuum,  (n«)  is  the  mean  free  electron  density  in  the  ionosphere,  and  re  is  the 
classical  electron  radius.  Under  the  small  angle  scattering  assumption,  it  can  be  shown  that  U obeys  the 
parabolic  wave  equation  [Tatarskii  1971]: 
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where  ne(r,t)  is  the  free  electron  density  which  is  a function  of  space  and  time. 
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Figure  1.  Propagation  geometry. 
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The  statistical  nature  of  n«(r,t)  requires  a statistical  approach  to  solving  the  parabolic  wave 
equation.  Thus  a differential  equation  for  the  mutual  coherence  function, 

I‘(Ar4co.At)  - <U(r1,(01,t1)U*(ir2,co2,t2))  (4) 

is  derived  from  the  parabolic  wave  equation  assuming  small  angle  scattering,  wide-sense  statistically 
stationary  processes,  and  a model  for  the  temporal  fluctuations  of  n*.  The  details  of  this  derivation  may 
be  found  in  Tatarskii  [1971]  or  Dana  [1986]. 

The  solution  of  the  differential  equation  for  the  mutual  coherence  function  depends  on  structure 
function  A(p)  [Sreenivasiah,  Ishimaru,  and  Hong  1976]: 


A(p) 


1 


exp  (iKi-p)  Se(K±,K2=0) 


d2K± 

(2n)2 


(5) 


where  p is  a position  vector  in  a plane  normal  to  the  line-of-sight  with  components  x and  y,  and  Se(K) 
is  the  three-dimensional  power  spectral  density  of  the  electron  density  fluctuations  in  the  ionosphere. 
Hie  structure  function  is  then  the  correlation  of  the  electron  density  fluctuations  in  a plane  normal  to 
the  propagation  direction.  It  is  these  fluctuations  that  cause  diffraction  and  scintillation  in  the  received 
signal.  A power-law  form  for  Se  is  usually  assumed  [Wittwer  1979]  which  results  in  a Bessel  function 
form  for  the  structure  function.  Because  the  structure  function  occurs  in  the  exponent  of  the  solution 
for  the  mutual  coherence  function,  the  quadratic  phase  approximation  is  made  to  produce  a 
mathematically  tractable  result.  This  quadratic  phase  approximation  simplifies  the  expression  for  A(p) 
to 


A(p)  = Aq  1 


A2 


(6) 


The  constants  Aq  and  A2  are  determined  by  functional  form  of  Se  and  by  the  phase  variance  imparted 
on  the  wave  in  the  ionosphere.  The  quantities  Lx  and  Ly  are  electron  density  fluctuation  scale  sizes  in 
the  plane  normal  to  the  line-of-sight. 


Although  much  information  is  available  on  the  spatial  fluctuations  of  ne  (see,  for  example. 
Wittwer  1979),  there  is  little  information  on  small  time-scale  fluctuations  of  ne.  Thus  S6  does  not,  in 
general,  contain  Doppler  frequency  information,  and  a model  is  required  to  describe  the  temporal 
fluctuations.  A general  model  for  temporal  dependence  of  the  quadratic  form  of  the  structure  function 
is 


A(p,t)  = Ao  1 


- A 2 


£ 
T 2 
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Cxt  xt 
LXT0 
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LyT0 


(7) 


The  quantity  T0  can  be  thought  of  as  a “decorrelation  time”  of  the  electron  density  fluctuations  in  the 
ionosphere.  The  degree  of  space-time  correlation  is  determined  by  the  coefficients  Cxt  and  Cyt  . 

This  form  for  the  structure  function  reduces  to  the  frozen-in  model  if,  for  example,  CXt  is  set  to 
unity  and  Cyt  is  set  to  zero.  Then  the  x dependence  can  be  written  as  (x-vt)2/Lx  where  the  velocity  is 
equal  to  LJT0,  The  turbulent  model  corresponds  to  the  case  where  both  Cxt  and  Cyt  are  zero. 

Given  a quadratic  form  for  the  structure  function,  a closed-form  solution  can  be  obtained  for  the 
mutual  coherence  function  [Knepp  1983,  Dana  1986].  This  solution  reduces  to  a mathematically 
tractable  form  when  it  is  assumed  that  the  scattering  occurs  in  a thin  layer.  With  these  assumptions,  the 
mutual  coherence  function  for  the  general  model  is 
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where  041  is  the  standard  deviation  of  the  phase  fluctuations  imparted  on  the  wave  in  the  ionosphere,  (do 
is  the  angular  frequency  of  the  RF  carrier,  and 


21}  Ik 


Ax  = U4x  + l4y. 


r 2ft*  ik 

“ U4x  + *y.  * 


The  decorrelation  distances  lx  and  4y  are  given  by  the  expressions 


(zt  + Zf)Lx 

WaT 


(Zt  + Zf)Ly 

zi^VaT 


the  decorrelation  time  is  given  by  the  expression 


T°  aoVXT  ’ 

and  the  coherence  bandwidth  <dcoh  is  given  by  the  expression 


AyCoftLj;  zt+  zr 

wcoh  - ~Ji , 


Wcoh  " 2c<4a2  *t*r  • 

GENERALIZED  POWER  SPECTRAL  DENSITY 

The  generalized  power  spectral  density  S(Ki,t,»od)  of  the  signal  incident  on  the  plane  of  the 
receiver  is  the  Fourier  transform  of  the  mutual  coherence  function: 


S(Ki,t,(0D)  = I d2p  J ^ J dt  r(p,(d,t)  exp  [-i(Kj.-p-(OT-(ODt)J 


where  angle-of-arrival  Ki  is  the  Fourier  transform  pair  of  position  p in  the  x-y  plane,  delay  x is  the 
Fourier  transform  pair  of  relative  carrier  frequency  co,  and  Doppler  frequency  <dd  is  the  Fourier 
transform  pair  of  relative  time  t,  The  quantity  S(Kx.t,(OD)(d2Ki/4n2)dT(d(OD/2n)  is  equal  to  the  mean 
signal  power  arriving  with  angles-of-arrival  in  the  interval  KJ4n2  to  (Ki+d2Ki)/4«2,  with  delays 
relative  to  a nominal  propagation  time  in  the  interval  x to  x+dx,  and  with  Doppler  frequencies  in  the 
interval  (0D/2n  to  (<dd+  dO)D)/2n.  The  delay  dependence  of  the  GPSD  is  a consequence  of  die  fact  that 
some  of  the  signal  energy  takes  a dog-leg  pfeth  through  the  ionosphere  from  the  transmitter  to  the 
receiver  and  arrives  later  than  the  signal  energy  that  propagates  straight  through  the  ionosphere. 
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In  general,  the  GPSD  can  be  written  as  the  product  of  a Doppler  spectrum  Sd(®d)  and  an  angle- 
delay  spectrum  SK-t(Kj.,T): 

* S(Kx,X,(0D)  = Sd(£0D)  SKTfw\±,T)  . (14) 

After  performing  the  integrals  indicated  in  Equation  13,  the  Doppler  spectrum  for  the  general  model  is 


Sd(»d)  = 


j/rrtjL 


Vl  - Cx2,  - C2 


exp 


(XqCI)D  - CxlKn^X  * CylKyty)2! 

4(1  - Cx2t  - Cyt)  J 


(15) 


In  terms  of  the  components  of  Kj.  (Kx  and  Ky),  the  angle-delay  part  of  the  GPSD  is 

W1?  r k y i y 1 

SKt(Kx,Ky,t)  = |jJ  2 My  OUOcoh  exp  [-  J 
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where  the  delay  parameter  a is  defmed  to  be 


a 


«o 


O^Otcoh 


(16) 


(17) 


The  value  of  a is  quite  large  under  strong  scattering  conditions.  The  components  of  Kjl  are  related  to 
the  scattering  angles  6X  and  0y  about  the  x and  y axis  respectively  by  the  relations 


2 n sin(Ox) 


2n  sln(6y) 


(18) 


It  should  be  noted  that  the  range  of  delay  in  Equation  16  is  from  -« to  +«•.  However,  this  delay 
is  relative  to  some  nominal  propagation  time,  and  the  value  of  the  GPSD  rapidly  approaches  zero  for 
decreasing  negative  values  of  <oCoh't-  In  the  limit  that  a is  infinity,  the  GPSD  is  non-zero  only  for 
positive  values  of  delay.  Thus  Equation  16  presents  no  real  problem  with  causality. 


It  is  interesting  to  examine  the  limits  of  the  general  model  Doppler  spectrum  in  order  to  show 
that  this  model  does  indeed  encompass  both  the  frozen-in  and  fully  turbulent  models.  These  limits  for 
the  Doppler  spectra  are: 


Limit 

Cxt-*1  Sd(gjd)  s 2n  to  -IMx) 

Cy(  — * 0 


“o  SD((oD)  * V^oexpf-^1  . 
Cyt-0  L 4 J 


(Frozen-in  Model) 


(Turbulent  Model) 


(19) 


(20) 


For  the  frozen-in  model,  the  delta-function  relationship  between  Doppler  frequency  and  Kx  is  what  is 
obtained  by  assuming  that  the  random  diffraction  pattern  of  the  signal  is  “frozen"  and  drifts  in  the  x 
direction  past  the  receiver.  For  the  turbulent  model,  the  Doppler  spectrum  is  independent  of  Ki  so  the 
temporal  and  spatial  variations  in  the  received  signal  are  also  independent. 


When  the  two  decorrelation  distances  are  equal  (lx  = ly  = l0),  the  scattering  is  isotropic  about 
the  line-of-sight,  and  the  GPSD  takes  a somewhat  simpler  form.  A three-dimensional  plot  of  the 
isotropic  one-dimension  angle-delay  part  of  the  GPSD  (Equation  16  integrated  over  one  component  of 
Ki)  is  shown  in  Figure  2.  This  plot  shows  the  mean  received  power  density  as  a function  of 
normalized  angle  Kl0  and  normalized  delay  (Ocoh?.  The  vertical  axis  is  linear  with  arbitrary  units. 


Figure  2.  Angle-delay  generalized  power  spectral  density. 


It  can  be  seen  that  the  power  arriving  at  large  angles  is  also  the  power  arriving  at  long  delays. 
The  power  arriving  at  long  delays  thus  has  higher  spatial  frequency  components  than  power  arriving  at 

short  delays.  When  there  is  strong  space-time  correlation  (i.e.  when  “v/cj,  + Cyt  is  approximately 
equal  to  unity)  these  higher  spatial  frequency  components  correspond  to  higher  Doppler  frequency 
components.  The  signal  arriving  at  long  delays  then  varies  more  rapidly  in  time  than  the  signal 
arriving  at  short  delays. 

Another  view  of  the  GPSD  can  be  obtained  by  considering  the  delay-Doppler  scattering  function 
(Equation  l'l  integrated  over  Ki).  For  the  turbulent  model,  the  scattering  function  is  separable  into  a 
function  of  Doppler  frequency  times  a function  of  delay.  This  is  rot  the  case  for  the  frozen-in  model. 

A comparison  of  the  scattering  functions  for  the  frozen-in  and  turbulent  models  is  shown  in 
Figure  3.  The  frozen-in  scattering  function  is  just  a reproduction  of  Figure  2 with  normalized  angle 
Kl0  replaced  with  normalized  Doppler  frequency  t0<dd.  This  is  a consequence  of  the  delta-function 
relationship  between  angle  and  Doppler  frequency  for  the  frozen-in  model.  For  this  model  the  signal  at 
long  delays  has  correspondingly  large  Doppler  shifts,  and  a wing-like  structure  is  seen  in  the  scattering 
function.  The  turbulent  model  scattering  function  does  not  exhibit  these  Doppler  wings  because  the 
Doppler  spectrum  is  the  same  at  all  delays.  Both  functions  have  exactly  the  same  power  density  at  each 
delay.  The  difference  in  appearance  of  the  figures  is  due  to  the  fact  that  the  turbulent  model  signal  at 
long  delay  is  more  spread  out  in  Doppler  frequency  and  is  therefore  less  obvious. 
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mOZEN-IN  MODEL  TURBULENT  MODEL 


Figure  3.  Comparison  of  scattering  functions  for  the  frozen-in  and 
turbulent  models. 


A progression  of  scattering  functions  for  the  general  model  is  shown  in  Figure  4.  The  space- 
time  correlation  coefficient  C*t  varies  from  0.99  for  the  scattering  function  in  the  upper  left  to  0.7  for 
the  scattering  function  in  the  lower  right.  The  scattering  function  for  Cxt  equal  to  0.99  is  essentially 
identical  to  that  for  the  frozen-in  model,  and  the  scattering  function  for  Cxt  equal  to  0.7  is  essentially 
identical  to  that  for  the  turbulent  model.  For  intermediate  values  of  Cxt,  the  scattering  functions  still 
exhibit  Doppler  wings  but  the  wings  have  broader  Doppler  spectra  as  Cxt  decreases. 


Figure  4.  Scattering  functions  for  the  general  model. 
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MATCHED  FILTER  EXAMPLE 


Signal  scintillation  can  be  either  flat  fading  (nonselective)  or  frequency  selective,  depending  on 
the  signal  bandwidth  and  the  degree  ionospheric  disturbance.  If  all  frequency  components  of  the 
received  signal  vary  essentially  identically  with  time,  the  propagation  channel  is  referred  to  as 
nonselective  or  flat  fading.  When  the  scintillations  exhibit  statistical  decorrelation  at  different 
frequencies  within  the  signal  bandwidth,  the  channel  is  referred  to  as  frequency  selective.  The  level  of 
frequency  selective  scintillation  disturbance  is  measured  by  the  frequency  selective  bandwidth,  f0.  This 
parameter  is  defined  as 


f _ I_  _ (Ocoh 

0 ~ 2itOx  - 2n 


(21) 


where  ax  is  the  standard  deviation  of  the  time-of-arrival  jitter. 

Techniques  have  been  developed  to  generate  representative  signal  structures  using  statistical 
channel  simulators  [Wittwer  1980  and  Dana  1986].  The  latter  reference  describes  techniques  for 
generating  realizations  of  the  channel  impulse  response  function,  and  describes  how  these  realizations 
can  he  used  to  generate  the  output  of  a matched  filter.  A report  is  currently  in  preparation  [Dana  1990] 
that  extends  these  techniques  to  inlcude  the  general  model. 

A comparison  of  the  matched-filter  output  amplitude  generated  with  the  frozen-in,  general,  and 
turbulent  models  is  shown  in  Figure  5 for  one  level  of  frequency  selective  propagation  disturbances, 
characterized  by  the  ratio  of  the  frequency  selective  bandwidth  f0  to  the  transmitted  chip  rate  Rc  (f(/Rc 
a 0.1  for  these  examples).  Each  frame  provides  a three  dimensional  picture  of  the  matched  filter 
output  for  a single  transmitted  pulse  as  a function  of  time  delay  (horizontal  scale)  and  time  (scale 
directed  into  the  figure).  The  total  duration  of  each  of  the  frames  is  10  decorrelation  times. 

The  top  frame  shows  an  example  generated  using  the  frozen-in  model  (Cxt  = 1 and  Cyt  = 0).  An 
effect  of  the  frozen-in  model  that  is  evident  in  the  top  frame  is  that  the  signal  arriving  at  long  delays 
varies  more  rapidly  in  time  than  the  signal  arriving  at  shorter  delays.  The  middle  frame  is  a general 
model  realization  (Cxt  = 0.9  and  Cyt  = 0),  and  the  bottom  frame  is  for  the  turbulent  model  (Cxt  = Cyt  ~ 
0).  The  difference  between  the  top  and  bottom  frames  is  that  the  turbulent  model  amplitude  has  the 
same  fading  rate  at  all  delays.  It  can  be  seen  that  the  general  model  realization  falls  somewhere 
between  these  two  limiting  cases. 
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DELAY  (CHIPS) 

Figure  5.  Matched  filter  output  amplitude  time  histories  for  frequency 

selective  channels. 
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